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Abstract 

The quantum discord is used as measure of quantum correlations for two families of multipartite 
coherent states. The first family interpolates between generalized GHZ states and generalized Werner 
states. The second one is an interpolation between generalized GHZ and the ground state of the 
multipartite quantum system. Two inequivalent ways to split the system in a pair of qubits are 
introduced. The explicit expressions of quantum quantum discord in multipartite coherent states are 
derived. Its evaluation uses the Koashi- Winter relation in optimizing the conditional entropy. The 
temporal evolution of quantum correlations (quantum discord and entanglement) is also discussed. 
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1 Introduction 



The characterization of quantum correlations in a bipartite or multipartite quantum system has been 
intensively investigated in the context of quantum information science [TJ [2] . Different methods for 
quantifying these correlations were reported in many works (for a recent review see [3]). Entanglement, 
one of the prominent of quantum correlations, constitutes a typical resource to manage information 
in several ways [H [5]. In this respect, a considerable attention has been paid to develop a quantita- 
tive theory of entanglement (concurrence, entanglement of formation and linear entropy [HI El [H])- 
Entanglement was for a long time viewed as synonymous with quantum correlation. However, some 
recent studies showed that entanglement is only a special kind of quantum correlations. Indeed, 
unentangled quantum states can also possess quantum correlations which play a relevant role in im- 
proving quantum communication and information protocols better than their classical counter-parts 
|10|-|15]. Therefore, as the non-classicality of correlations present in quantum states is not due solely 
to the presence of entanglement, there was a need of a measure to characterize and quantify the non- 
classicality or quantumness of correlations which goes beyond entanglement. The first attempt in this 
direction was made in the works [16] and [T7J where the authors concluded that when entanglement is 
subtracted from total quantum correlation, there remain correlations that are not entirely classical of 
origin. Nowadays, it is commonly accepted that the most promising candidate to measure quantum 
correlations is the so-called quantum discord. It has attracted considerable attention and continues to 
be intensively investigated in many contexts such as quantum decoherence, quantum computation and 
phase transition as well as other related fields (a detailed list of references can be found in [3]). Many 
efforts were deployed to get a closed analytical expression of quantum discord. Only some partial 
results were derived for some special two-qubit systems [16|. I1T |. 1X8] , I19j . The extension of the notion of 
the discord to continuous variable systems and the discussion of its properties were considered in [20J . 
It must be emphasized that there is now an intense recent research activity to demonstrate experimen- 
tally the advantages of the use of quantum discord without entanglement in quantum protocols. In a 
recent work [21], the authors experimentally show, using a variety of polarization-correlated photon 
pairs, non-zero quantum discord is the optimal resource for remote state preparation. More interest- 
ingly, it is found that unentangled state with non-zero quantum discord provides better performance 
than entangled state by achieving a higher fidelity (which is directly related to quantum discord) in 
remote state preparation. 

In this paper we shall be mainly concerned with the pairwise quantum correlations, especially 
quantum discord, present in multipartite coherent states. It must be noticed that the entangled co- 
herent states have received a special attention in the last two decades (for a recent review, see [22]). 
Therefore, paralleling the treatment of entanglement in coherent states systems, the main of this work 
is to investigate the pairwise quantum discord present in multipartite nonorthogonal states. In fact, en- 
tangled coherent states, which are typical examples of entangled nonorthogonal states, have attracted 



2 



much attention in the literature. One may quote for instance entangled Glauber coherent states [23J, 
SU{2) and SU(1, 1) entangled coherent states [24] , The term entangled coherent state was introduced 
by Sanders in a study concerning production of entangled coherent states by using a nonlinear Mach- 
Zehnder interferometer [23\ I25j . Entangled coherent states were initially treated as bimodal states 
but later generalized to superpositions of multimode coherent states [261 EH 125] - Generalizations to 
multimode systems allow the intricacies of multipartite entanglement to become manifest in entangled 
coherent states, as for instance in GHZ (Greenberger-Horne-Zeilinger), W (Werner) states [291130] and 
entangled coherent state versions of cluster states [3111321 1315] . In this sense, the characterization of the 
quantum discord in nonorthogonal multipartite states constitutes an important issue which deserves 
the same degree of attention as entanglement. 

This paper is organized as follows. In the first section, we introduce a special instance of superpo- 
sitions of multipartite coherent states. They involve Weyl-Heisenberg, SU{2) and SU(1, 1) coherent 
states labeled by a single complex parameter z. We shall focus on balanced superpositions which 
are symmetric or antisymmetric under the parity transformation z — > —z. The symmetric super- 
position is interpolating between the ground state of the multipartite system and GHZ state. The 
antisymmetric superposition provides a continuous interpolation between the Werner state and GHZ 
state. To study the bipartite quantum correlations present in the system, two different qubit mapping 
are considered. Section 3 is devoted to the derivation of the explicit expression of the quantum dis- 
cord. The optimization over all the measurement is performed by combining the purification method 
and the Koashi- Winter relation . Finally, in the last section we discuss the dynamical evolution of 
entanglement and quantum discord under a dephasing channel. Concluding remarks close this paper. 

2 Multipartite coherent states and qubit mapping 

For several quantum systems, the coherent states can be obtained by exploiting the structure relations 
of the relevant dynamical group structure G. Therefore, any quantum state in the corresponding 
Hilbert space can be expanded as a sum of the coherent states associated with the group G. In 
this respect, a widespread interest was devoted to the theoretical as well experimental studies of 
superpositions of coherent states. Evidences of such superpositions first appeared in a study of a 
certain type of nonlinear Hamiltonian evolution in [34} 135]. Also, a detailed analysis concerning the 
manifestation of superpositions of coherent states was reported in [361 137] (see also [38]). However, it 
is important to stress that superpositions of coherent states are experimentally difficult to produce, 
and fundamentally this could be due to extreme sensitivity to environmental decoherence. Some 
experimental efforts to create superpositions of coherent states were reported in [39]. 
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2.1 Basic of coherent states 

A system of coherent states is defined to be a set £ T>} of quantum states, in some Hilbert 

space H, parameterized by the variable z belonging to the set V such that: z — > \z) is continuous 
and the system is over-complete; i.e. 

J '\z){z\ dfi(z) = I H . (1) 

The continuity and the resolution of the unity ([I]) form the minimal set of requirements to define a 
system of coherent states. They are not orthogonal to each other with respect to the positive measure 
in (H|). It should be noticed that the explicit form of the coherent state \z) in not required for our 
purpose. However, to illustrate some interesting limiting situations and for the sake of completeness, 
we consider the familiar sets of Perelomov coherent states [40] associated with Weyl-Heisenberg, SU(2) 
and SU(1, 1) groups (the extension to other groups and other coherent states definitions is straight- 
forward). A simple way to deal with these symmetries in a unified scheme can be achieved by using 
the so-called generalized Weyl-Heisenbeg algebra A (see |41j and references therein) spanned by an 
annihilation operator a~ , a creation operator a + and a number operator N ^ a + a~ satisfying the 
structure relations 

[a-,a + ] = G(N), [N,a~} = -a~, [N,a + ] = +a + , (2) 

with 

o + = (a~) t , N = N*, G(N) = F(N + 1) - F(N), (3) 

where the F structure function characterizes the deviation from the usual harmonic oscillator. A 
representation of the algebra A which extends that of Weyl-Heisenberg, is defined through the actions 

ar\n) = y/F(n)\n-l), a~|0) = 0, (4) 
a + \n) = ^/F(n + l)\n + 1), N\n) = n\n) (5) 

on the Hilbert space spanned by the eigenvectors of N . This formally defines a representation of A. 
Note that a + a~ = F(N), a relation that generalizes N = a + a~ for the harmonic oscillator. The 
operator F(N) can be considered as the Hamiltonian for a quantum system. This provides a physical 
significance to the structure function F. The representation (Fock-Hilbert) space % is generated by 
the orthonormal set {|n) : n ranging}. The dimension of the representation afforded by (HD and (0) 
is controlled by the positiveness of the eigenvalues F{n) of the operator F{N). The usual Weyl- 
Heisenberg algebra is recovered for 

F(N) = N. 

The algebra A coincides with su(l, 1) Lie algebra for 

F(N) = N(2k -1 + N) 
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where the number k, which acquires discrete values k = |,2, characterizes the SU (1,1) 
discrete-series representations. The algebra A reproduces the spin algebra su(2) when 

F(N) = N(2j + 1-N), 

and the creation and annihilation generators coincides with the usual raising and lowering operators 
acting in the spin representation of dimension 2j + 1. 

The Perelomov coherent states for the algebra A can be constructed via a simple strategy based 
on the use of a Fock-Bargmann space associated with A [41J. This is sketched briefly as follows. Let 
us look for states in the form 

\z) =^2a n z n \n), a„£l, z £ C, (6) 
n 

where the sum on n is finite or infinite according to as A admits a finite- or infinite-dimensional 
representation. The a n coefficients can then be determined from the correspondence rules 

|n) — ► a n z n , aT — > (7) 

applied to relations (J1J) and ([5]). This analytical realization leads to the following recurrence relation 



na n = y/F(n)a n -i, (8) 

which can be iterated to give 



n! 

(by taking ao = 1). In Eq. 0, the generalized factorials are defined by 

F(0)! = 1, F(n)\ = F(1)F(2) . . . F(n). 
This yields the following normalized coherent states 



(9) 



\ z )=N(\z\)^V^z n \n) (10) 



where the normalization factor is 

-2 



N(\z\) 



E 



(n!) 2 11 ' 



subject to convergence. They satisfy \z) = exp(za + )|0) and are thus coherent states in the Perelomov 
sense. For the usual harmonic oscillator (F(n) = n), the equation (|1Q[) reads 



, v _Mi \ - z 

\z) = e 2 



n=0 v'm 
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The overlap between two Glauber coherent states is 

( Zl \z 2 ) = exp(--(|2i| 2 + \z 2 \ 2 - 2z lZ2 )). (12) 

The standard set of the SU(2) coherent state is obtained from (jlOp for F(n) = n(2j + 1 — n). They 
are given by 



2j 



nl/2 



n\(2j — n)\ 



z n \n). (13) 



n=0 

The parameter z can acquire any complex value. The SU{2) coherent states are normalized but they 
are not orthogonal to each other: 

( Zl \z 2 ) = (1 + M V(l + N V'U + ziz 2 ) 2j . (14) 
Finally, for F(n) = n(2k — 1 + N), one gets the SU(1, 1) coherent states 

r (2/fc-l + n)l" 



\z) = (i-\ z \rE 



nl/2 



n\(2k - 1)! 



z n \n). (15) 



n=0 

where the complex variable z belongs to the unit disc {|z| £ C, |z| < 1}. The kernel of two SU(1, 1) 
coherent states reads 

( Zl \z 2 ) = (1 - M 2 ) fc (l - \z 2 \ 2 ) k (l - z lZ2 )- 2k . (16) 
2.2 Superpositions of multipartite coherent states 

The over-completeness relation makes possible the expansion of an arbitrary state of the Hilbert space 
"H in terms of the coherent states of the system under consideration. It follows that when one considers 
a collection of n noninteracting identical particles, the whole Hilbert space is a tensor product and 
any multipartite coherent state can be written as a superposition of tensorial products of the form 
\ z \) <8> \z 2 ) ■ ■ ■ <S> \z n ) = \zi, z 2 , ■ ■ ■ z n ). In fact, the resolution of identity allows us to expand any state 
l^) in the space H ® T~L <g> • • • H as 

dfi(zi) d[i(z 2 ) ■ ■ ■ dfj,(z n )\z 1 ,z 2 , ■ ■ ■ z n )(zi, z 2 , • • • z n \ip) (17) 

reflecting that any multipartite state can be viewed as a superposition of the coherent states \zi, z 2 , • ■ ■ z n ). 
The multipartite state (|17p can be reduced to a sum if the function ip{z\, z 2 , ■ ■ ■ z n ) = (z\,z 2 , ■ ■ ■ z n \ip) 
can be expressed as a sum of delta functions. Indeed, setting 

ip{zi,z 2 , ■■■z n ) = {5(z - zi) + e imn 5(z + z 1 ))5(z 1 - z 2 )5(z 2 - z 3 ) ■ ■ ■ 5(z n -i - z n ), 

one gets the following balanced or equally weighted superposition of multipartite coherent states 

|^) = \z,m,n) = M{\z) ® \z) ® ■ ■ ■ ® \z) + e im7T \ - z) ® \ - z) ® ■ ■ ■ ® | - z)) (18) 

where m G Z. The normalization factor N is 



M = [2 + 2p n cos 



mir\ 



-1/2 



where p denotes the overlap between the states \z) and | — z). It is given 

p=(z\- z). 

It is real as it can be verified from (|10p . For Weyl-Heisenberg, SU{2) and SU (1,1) coherent states, 
the quantity p is obtainable from the expressions (|12p . (|14p and (|16p . respectively. Two interesting 
limits of superpositions of the form (|18p arise when p — > and p — > 1. 

We first consider the asymptotic limit p — > 0. In this limit the two states \z) and | — z) approach 
orthogonality, and an orthogonal basis can be constructed such that |0) = \z) and |1) = | — z). Thus, 
the state \z,m, n) approaches a multipartite state of GHZ type 

\z,m,n) ~ |GHZ) n = -^(|0) ® |0) ® • • • ® |0) + e im ^|l) ® |1) ® ■ ■ ■ ® |1}). (19) 
v2 

In the situation where p — > 1 (or z — >■ ), one should distinguish separately the cases m = (mod 2) 
and m = 1 (mod 2). For m even, the multipartite superposition (|18p reduces to ground state 

|0,0 (mod 2),n) ~ |0) ® |0) ® ■ ■ ■ ® |0), (20) 

and for m odd, the state \z, 1 (mod 2), re) reduces to a multipartite state of Werner type [32] 

|0,1 (mod2),n) ~ |W>„ = -=(|1)® |0)®- • -® |0) + |0)® . .® |0) + - • - + |0)® |0)®- • . (21) 

y/n 

Here |n) (n = 0, 1) denote the Fock-Hilbert states. 

It follows that the states \z, m = (mod 2), n) interpolate between states of GHZ type (p — > 0) and 
the separable state |0) ® |0) ® • • • ® |0) (p — > 1). In other hand, the states \z, m = 1 (mod 2), re) may 
be viewed as interpolating between states of GHZ type (p — > 0) and states of Werner type (p — > 1). 

To close this subsection, it is important to emphasize that the main concern in investigating the 
properties of superpositions of coherent states is how to produce such states. Their experimental 
production is fundamentally difficult to achieve. This is especially due to extreme sensitivity to 
environmental decoherence. Experimental efforts to create superpositions of coherent states with the 
present day technology, are encouraging. Indeed, superpositions of weak coherent states with opposite 
phase, resembling to a small "Schrodinger's cat" state (or "Schrodinger's kitten"), were produced by 
photon subtraction from squeezed vacuum [33] . Also, the experimental generation of arbitrarily large 
squeezed Schrodinger cat states, using homodyne detection and photon number states (two photons) as 
resources was reported in [43]. Very recently, creation of coherent state superpositions, by subtracting 
up to three photons from a pulse of squeezed vacuum light, is reported in [33]. The mean photon 
number of such coherent states produced by three-photon subtraction is of 2:75. The production of 
cat states especially ones of high amplitude or mean number of photons remains an experimental 
challenge. Considering the fast technical progress and the increasing number of groups working in this 
field, we expect that the generation of cat states (and Bell states) is a goal that is achievable in the 
near future. 
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2.3 Qubit mapping 
2.3.1 Pure states 



To study the bipartite quantum correlations present in (|18p . the whole system can be partitioned in 
two different ways. We first consider bipartite splitting of the multipartite system, i.e., splitting the 
entire system into two subsystems, one subsystem containing any k (1 < k < n — 1) particles and the 
other containing the remaining n—k particles. Accordingly, one writes the state \z,m,n) as 

\z,m,n) = M(\z) k ® |z) n _ fc + e imn \ -z) k ®\- z) n - k ) (22) 

where 

| ± z)i = | ± z)i ® | ± z)2 <8) • • • <8) | ± z);, I = k,n — k. 

The multi-particle state |z,m,n) can be expressed as a state of two logical qubits. For this end, we 
introduce, for the first subsystem, the orthogonal basis {|0)fc, |l)fc} defined as 

|Q) fc = M^L^ |iv -!*>*" I "*>* (23) 
V2(l+P*) V2(l"P fc ) 

Similarly, we introduce, for the second subsystem, the orthogonal basis {|0) n _fc, |l)n-fc} given by 

\ n \ \ z )n—k + I ~~ z )n~k \-.\ \ z )n— k ~ I — Z) n —k . . 

U/n-fc = > ; = — M-M-Jfc = / = = — ■ (24J 

V2(l+p"~ fc ) ^2{l-p n - k ) 
Reporting the equations (|23p and (|24p in (|22p . one has the explicit form of the pure state \z,m,n) in 
the basis {|0) fc ® \0) n -k, \0)k ® |l)n-Jfc, ® |0) n _ fc , ® |l) n -fc}- It is given by 

\z,m,n) = Ca 

(25) 

o=0,l/3=0,l 

where the coefficients C a ,f3 are 

Q),o = M{1 + e im > fe a n _ fe , C ,i = M{1 - e im7 > fe 6 n _ fe 

Ci )0 = AA(1 - e^)a n _ fc 6 fc , Ci,i = N{1 + e^)b k b n . k . 
in term of the quantities 



ai = V 1 ~2 ^ ' bl = y <l for I = k,n- k 

involving the overlap p between two coherent states of equal amplitude and opposite phase. 

2.3.2 Mixed states 

The second partition can be realized by considering the bipartite reduced density matrix p k \ which is 
obtained by tracing out all other subsystems except ones labeled by the indices k and I. There are 
n(n — l)/2 different density matrices p k \. It is simple to see that all the reduced density matrices p k \ 
are identical. Therefore, it is sufficient to consider p\2 and to generalize from this case. Then, by 
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tracing out systems 3, 4, . . . , n in the state \z, m, n), we obtain the reduced density matrix describing 
the particles or modes 1 and 2 as 



pu = r Tt 3A: ,„ :n (\z,m,n)(z,m,n\) 

= Af 2 (\z, z)(z, z\ + | - z, -z)(-z, -z\ + e 1 



z, —z)(z, z\ + e 



\z,z)(-z,-z\) (26) 



with q = p n 2 . To study the correlations of the system described by the density matrix pi2, we convert 
it into a two-qubit system. Thus, we choose an orthogonal pair {|0), |1)} as 



\z) = a\0)+b\l) 



z) = o|0) - 6|1) 



(27) 



where 



1 +p 



1 — p 



The logical qubits |0) and |1) given by 

1 







v / 2T2^ 



(l*> 



*)) 



i 



(l*> 



z)). 



(28) 



coincide with even and odd coherent states, respectively. Substituting the equation (|27jl into (1261) . we 
obtain the density matrix 



Pi2 = A/" 2 



2a 4 (l+gcosm7r) 2a 2 b 2 (l+q cos run) ' 

2a 2 fe 2 (l— gcosm7r) 2a 2 fe 2 (l— q cos m7r) 

2a 2 b 2 (l— gcosm-7r) 2a 2 b 2 (l— q cos m7r) 

2a 2 6 2 (l+qcosm7r) 2b 4 (l+g cos mn) 



(29) 



in the basis {|00), | OX) , |10), | XX) }-. It is remarkable that the obtained density belongs to the set of 
the so-called X states. It can be also written, in the Bloch representation, as 

3 3 
0=0/3=0 

where the correlation matrix R is given by 

/ 1 2M 2 p(l + qcosrrnr) \ 

2AA 2 (l-p 2 ) 

-2A/" 2 (1 -p 2 )qcosrmr 

V 27V 2 p(l + gcos?n7r) 2J\f 2 (p 2 + q cos mir) J 

and <r° and a 1 {% = 1,2,3) stand for identity and the usual Pauli matrices, respectively. 



R 



(31) 



3 Quantifying the quantum discord 

So far, quantum discord has been calculated explicitly only for a rather limited set of two-qubit 
quantum states and analytical expressions for more general quantum states are not known. In this 
section, we give another instance of states for which quantum discord can be evaluated explicitly. 
Indeed, we shall derive the explicit form of this kind of bipartite correlation in the state (|18p by 
making use of the qubit mapping corresponding to the partitioning schemes ([22]) and ([26]) discussed 
in the previous section. 
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3.1 Definitions 

The quantum discord is denned as the difference between total correlation and classical correlation. 
The total correlation is usually quantified by the mutual information I 

Hpab) = S( PA ) + S( PB ) - S{ PAB ), (32) 

where pab is the state of a bipartite quantum system composed of the subsystems A and B, the 
operator pa(b) = ^b(A){pab) is the reduced state of A(B) and S(p) is the von Neumann entropy of a 
quantum state p. Suppose that a positive operator valued measure (POVM) measurement is performed 
on subsystem A. The set of POVM elements is denoted by M = {M k } with M k ^ and M k = I. 
In this paper, we deal with two-qubit rank two states (with only two nonzero eigenvalues). In this 
case, the generalized positive operator valued measurement is not required. Indeed, it has been clearly 
shown in |45} HB] (see also the recent review [37] ) that for a bipartite mixed state of two qubits of rank 
two , the optimal measurement giving the quantum discord is a two element POVM. The elements of 
such POVM are orthogonal projectors. Thus, it is natural to consider only projective measurements 
for the subsystem A. The von Neumann measurement (from now on just a measurement) on the 
subsystem A project the system into a statistical ensemble {pB,k, PB,k} such that 

(M k ®I)p AB (M k ®I) 

PAB > 

PB,k 

where the measurement operation is written as [18] 



M k = UIl k tf (33) 

with U k = \k)(k\ (k = 0,1) is the projector for subsystem A along the computational base \k), 
U € SU(2) is a unitary operator and 

PB ,k = tr (M k <g> l)p AB (M k I) 

The amount of information acquired about particle B is then given by 

S{pb) ~ PB > k S (PB,k), 

k 

which depends on measurement Ai. This dependence can be removed by doing maximization over all 
the measurements, which gives rise to the definition of classical correlation 

C{pab) = max^ S(p B ) - J2k VB,k S{p Btk ) 

= S{p B ) - S min (34) 

where iSmm denotes the minimal value of the conditional entropy 

S = ^2 PB > k S (pB,k)- (35) 
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It follows that quantum discord is then given by 

D(pab) = I(pab) ~ C(pab) = S(p A ) + S min - S(pab)- (36) 

The main step in evaluating the quantum discord is the minimization of conditional entropy to get an 
explicit expression of the quantum discord in the multipartite system (|18p . 

3.2 Quantum discord for pure bipartite coherent states 



According to the first partitioning scheme (J22J), the entire n particles system is split into two com- 
ponents A, containing k particles, and B, containing n — k particles. The bipartite density state 
Pk,n-k = \z,m,n){z,m,n\ is pure and the conditional density is also a pure state. This implies that 
the quantum conditional entropy is zero. Then, the quantum discord for the pure state pab = Pk,n-k 
coincides with the von Neumann entropy of the subsystem A: 

D( Pk , n -k) = S( Pk ) (37) 

where pk is the reduced density of the subsystem A. In this scheme the quantum discord can be 
computed easily. It is given by 

D{p k ,n-k) = -A+ log 2 A+ - A_ log 2 A_ (38) 
in term of the eigenvalues of the reduced density matrix p k given by 



X± = l[l± x /1-Cl n _ k ) (39) 



where C k) n-k is the concurrence between the subsystems A and B: 

yj\ -pik^i _ p 2(n-k) 

Ck,n-k = 2|Co,oCi,i — Ci,oCo,i| = — — (40) 

1 + p n cos mir 



where we used the mapping (|25p . Note that the entanglement of formation given by 

E( Pt „- k ) = H(\ + f+f: t ™™ \ (41) 
\ 2 2 1 + p n cos run J 

is nothing but the von Neumann entropy of the subsystem A. Here H stands for the binary entropy 
H(x) = h(x) + h(l — x) with h (x) = —x log 2 x. Hence, we have the following closed relation between 
quantum discord and entanglement of formation 

D{p k , n _k) = E(p Kn _ k ). (42) 

This agrees with the fact that quantum discord and entanglement of formation are identical for pure 
states and amount to the same set of correlations. 

In the limiting case p — > 0, the state (|22p is of GHZ type having a maximal bipartite entanglement 
(Ck,n~k = 1) and we obtain 

D(p k ,n-k) = 1- 
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The situation becomes different when p — > 1. In this case, we have 

D{pk,n-k) = 

for symmetric pure states (i.e. m even) as expected (in this limit the state (|18p is a n tensorial product 
of the ground state |0)). In the limit p — > 1, the antisymmetric states (i.e. m odd) are of Werner type. 
The bipartite concurrence is 

Ck,n-k = ~\/ k ( n - *0> 

n 

and the corresponding pairwise quantum discord takes the form 

v TT (l l\n-2k\ 



2 2 n 

The quantum correlations (entanglement and quantum discord) in multipartite Werner states vanish 
as n becomes large. 

3.3 Quantum discord for mixed bipartite coherent states 
3.3.1 Mutual information entropy 

The subsystems A and B of the subsection 3.1 correspond here to the modes or particles 1 and 2 
described by the density operator p±2 (|29p . The non vanishing eigenvalues of the density matrix p±2 
are 

A± = 1 (l±p 2 )(l±qcos(m^)) ^ 

2 1 + p n COs(?7T,7r) 

and the joint entropy is 

S(p 12 ) = h(\+) + h(\-) = H(\+). (44) 

The eigenvalues of the marginal p\ = T1C2P12 are 

1 (1 ± p)(l ± pq cos(m7r)) 



A 



i,± 



2 1 + p n cos(m7r) 
and the marginal entropy reads 



S( Pl ) = h(X li+ ) + ^(Ai,_) = ff(Ai,+). (45) 



The eigenvalues of the marginal P2 = Tri/?i2 are 

1 (1 ± p)(l ± pq cos(m7r)) 



A 



2,ziz 2 1 + p n cos(m7r) 
and the corresponding entropy is given by 

S( P2 ) = h(X 2 ,+) + h(X 2 ,-) = H(X 2 ,+). (46) 

Note that for the mixed density under consideration the marginal densities pi and p2 are identical. 
The explicit form of the mutual information (|32p writes 

Ho ) - 2H ( 1 (1 + P) (1 + cos ( m7r )) ^ H ( 1 (1 + P 2 ) (1 + P"" 2 cos(^7t)) A 
W2J V 2 l+p n cos(m7r) / V 2 l+p n cos(mvr) / ' 1 J 
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3.3.2 Conditional entropy 

After computing the quantum mutual information, we next compute the classical correlation C{p\2) = 
C(pab)- This requires an explicit expression of the unitary operator U occurring in (|33j) . So, we 
parameterizes U as follows 

U = exp(r/cr + — fjcr-) exp(i(j)as) r) & C, (del. 

Using this parametrization, one can verify that the quantities defined by 

(<Ti) k = (k\U j aiU\k), i = 1,2,3 and k = 0, 1 

are given by 

(^)k = {-r—^, (<Ti)k = (-) k -—- Av2)k = i{-T-—- (48) 
1 + aa 1 + aa 1 + aa 

where a = —i-^L= tan y/fjrj. They can be also written as 

(<rs)k = (-) k cos6, (ax) k = (-) k sm9 cos tp , (a 2 )k = (-f sin^sin^ (49) 
where \e lLp = —irj. The conditional density operators ps,k = P2,k are 

P2 ' k P2,k { Rii(ax) k + iRnfah (1 - R03) + ORso " ^33)(^ 3 )fc J 



where the matrix elements R a p are giving by (|3ip and 

P2,k = + ^30 (<7 3 >fc). (51) 

It follows that the conditional entropy given by 

£ = V2,kS{p2,k) (52) 
fc=o,i 

can be cast in the following form 



S= P2,k n(± + ±y/l-4detp2, k \ 

fc=0,l ^ ' 



(53) 



Evidently, the quantities £>2,fc and det can be expressed in terms of the angular directions 6 and (p 
by making use of equations (f4U|) . ([50]) and (|5"T|) . Consequently, we can explicitly write S = S(6, (p) and 
perform the minimization over the azimuthal and polar angles 8 and tp. However, there exists another 
elegant method to perform such an optimization. This is presented in the following subsection. 
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3.3.3 Minimization of conditional entropy 



To minimize the conditional entropy (|53p , we shall use the purification method and the Koashi- Winter 
relation [48] (see also [49]). This relation establishes the connection between the classical correlation 
of a bipartite state pab = P12 an d the entanglement of formation of its complement pbc = P23- 
This connection will be clarified hereafter. The density pi2 (|29p is a two-qubit state of rank two. It 
decomposes as 

pi2 = X+\<f>+)((f>+\+X-\<l>-)((l>-\ (54) 

where the eigenvalues A+ and A- are given by (J3HJ) and the corresponding eigenstates \<t>+) and 
are given by 



1 



(l+p)|0,0) + (l-p)|l,l) 



1 

V2 



|0,1) + |1,0) 



(55) 



in the basis ([28]) . Attaching a qubit 3 to the two-qubit system 1 and 2, we write the purification of 
pu as 

D|l> (56) 
from which one 



such that the whole system 123 is described by the pure density state p±23 = 
has the bipartite densities p±2 = TV3/9123 and P23 = TV1P123. 

Suppose now that a von Neumann measurement {Mo, Mi} is performed on the subsystem 1 (here 
also we need positive operator valued measurement of rank one that is proportional the one dimensional 
projector). From the viewpoint of the whole system in the pure state \(j>), the measurement gives rise 
to an ensemble for P23 that we denote by 



where 



£23 = {Pk, \4>23,k)} 
\4>23,k) (<p23,k 



Tn 



Pit 



(M fc 01(8)1)1^)^1 



On the other hand, from the viewpoint of the state p±2, the von Neuman measurement on 1 gives 
rise to the ensemble for p2 that it was defined previously as £2 = {p2,k,P2,k} (note that pk = P2,k)- It 
is simple to check that the ensemble £2 can be induced from £23 by tracing out the qubit 3, namely 



P2,k 



Tr, 



|</>23,fc)(<?Wl 



We denote by E{\4>23,k)) the measure of entanglement for pure states. It is given by the von Neumann 
entropy of the reduced subsystem p 2 ,fc = Tr3(|023,fc)(023,fc|) 

E{\(t>23,k)) = S(p 2 , k ). 
It follows that the average of entanglement of formation over the ensemble £23 

^23 = Yl PkE{\<ha,k)) 

fc=0,l 
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coincides with the conditional entropy (|53|) . At this point, it is important to note that Koachi and 
Winter have pointed out that the minimum value of E23 is exactly the entanglement of formation of 
P23- Consequently, the minimal value of the conditional entropy coincides with the entanglement of 
formation of P23 ■ 

S m in = E(p 2 3), (57) 



which is easy to evaluate. Indeed, we get 



1 1 



S m m = E(p 23 ) = H(- + -Vl-|C(P2 3 )| 2 ) (58) 



where the concurrence of the density P23 is 

p 2 (l-p 2 )(l-p 2 "- 4 ) 



\C{P23)Y 




(1 + p n cos mir) 2 
Using the equation (|53p . one can verify the relation 

S min = S(9 = ^,0 = O) 

which reflects that the minimal value of conditional entropy is reached for 9 = ^ and = 0. Finally 
reporting (|44j) , (|45j) and (f58j) in the definition (f36l) , the explicit expression of quantum discord for the 
density p\2 is 

/ 1 (1 + p)(l + p"- 1 cos(m7r)) \ / 1 (1 + p 2 )(l + p n ~ 2 cos(mTr)) \ 
{Pl2) \2 l+p n cos(m7r) / V 2 l+p n cos(m7r) / 1 ' 

p 2 (l-p 2 )(l-p 2 "- 4 ) 
(1 + p n cos run) 2 

in term of the the overlap p. Clearly, this result can be also derived using the method developed 
(independently) in [19] and |50] (see also [51]) to obtain the quantum discord for some special instances 
of the so-called X states. However, one should stress that our method based on the purification 
trick combined with the Koashi- Winter relation reduce drastically the optimization process of the 
conditional entropy. It provides a direct relation between quantum discord and entanglement of 
formation reflecting how these correlations are distributed in a given tripartite pure system. It has 
been used in several papers dealing with quantum correlations as for instance in [52] where the authors 
emphasized its crucial role in exploring the distribution of entanglement in the deterministic quantum 
computation with one single pure qubit and a collection of an arbitrary number of mixed states. 

3.3.4 Some particular cases 

We start with the special case n = 2. The state ([26]) is pure and q = 1. The quantum discord (]59|) 
reads 

/ (1 +p)(l +pcosm7r) 

Diftia) = HI 2 r— 

\ 2(1 + p z cos mir) 
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It coincides with the quantum discord -D(/Jfc >n -fc) given by the equation (|4ip for n = 2 and k = 1. In 
particular, for symmetric states (m even), one obtains 

and D(p\2) = = 1 for antisymmetric states (m odd) as shown in the figure 1. 

Here also, it is interesting to consider the limiting cases p — > and p — > 1 as it was done previously 
for pure states. For n > 2, the quantum discord vanishes (|59p when p — > 0. It vanishes also when 
p — > 1 for m even. However, for m odd, the quantum discord (j59[) reduces to the special form 

and goes to zero for n large. 

To corroborate our analysis, we give in the figures 2 and 3 the behavior of the quantum discord 
(1591) . Figure 2 gives a plot of quantum discord versus the overlap p for the mixed state p\2 with m 
even (the symmetric case). As seen from the figure, after an initial increasing, the quantum discord 
decreases to vanish when p — > 1. The maximum of quantum discord depends on the number of particles 
contained in the system. It is remarkable that for n = 25 the maximum is larger than ones obtained 
for n = 4 and n = 5. Beside the numerical results reported in the figure 2, we have also studied 
the behavior of quantum discord for other values of n. This study shows that for 6 < n < 25 the 
maximum of quantum discord is greater than one reached for n = 5. Also, the behavior of quantum 
discord for n > 20 is very close to the case n = 25 presented in the figure 2. In figure 3, we give a 
plot of the quantum discord (I59p for m odd (the antisymmetric case) and different values of n. In this 
case the quantum discord increases as p increases for a small number of particles and the maximal 
value of quantum discord is reached for p — )■ 1. However, for a higher number of particles n (n = 25 
for instance), the maximum is reached for p < 1. In the limit p — > 1, we have a Werner state \W) n 
and the pairwise quantum discord given by (I60p decreases to vanish in the limit of a large number of 
particles (n — > oo). For p — > 0, corresponding to GHZ type states, the quantum discord is zero. 



D 




P 



FIG. 1: The pairwise quantum discord D versus the overlap p for n = 2. 
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FIG. 2: The pairwise quantum discord D versus the overlap p for symmetric states. 




FIG. 3: The pairwise quantum discord D versus the overlap p for anti-symmetric states. 

4 Dynamics of quantum correlations under dephasing channel 

The sudden disappearance of entanglement is one of the most intriguing features in quantum me- 
chanics. In fact, it has been observed that in a pair of entangled qubits, interacting with noisy 
environments, entanglement can disappear in a finite time [53 . This phenomenon, termed in the 
literature "entanglement sudden death", was experimentally confirmed |54j . 

In this section, we investigate the dynamics of bipartite quantum correlations (entanglement and quan- 
tum discord) of the multipartite coherent states \z, m, n) (n > 2). We focus on the second partitioning 
scheme (j26|) where mixed states emerge. The quantum discord and entanglement are two different 
correlations for mixed states (contrarily to the pure case) . We consider the evolution of a mixed bipar- 
tite system under a dephasing dissipative channel. In this order, we use the Kraus operator approach 
[3] which describes conveniently the dynamics of two qubits interacting independently with individual 
environments. The time evolution of the bipartite density p\2 = Pi2(0) (|29h can be written compactly 
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as 



p i 2(t)=Y,E lhV (t) p 12 (0) El„{t) 
where the so-called Kraus operators 

E*u (*) = E^t) ® E v {t) Y, E U E ^ = l - 



The operators E^ describe the one-qubit quantum channel effects. The non-zero Kraus operators for 
a dephasing channel are given by 

E = diag(l, Vl-7) E x = diag(0, ^7) 

with 7 = 1 — e~ Tt and T denoting the decay rate. It is easy to check that the density matrix (I29P 
evolves as 



2a 4 (l+<jcosm7r) 2(l--y)a 2 b 2 (l+q cos rmr) ' 

2a 2 i> 2 (l— qcosm7r) 2(l-7)a 2 b 2 (l-q cos m7r) 

2(l-7)a 2 fc 2 (l-gcosrmr) 2a 2 fe 2 (1-q cos rmr) 

2(l-7)a 2 6 2 (l+<?cosm7r) 2b 4 (l+<j cos m7r) 



(61) 



Using the prescription provided in the works |55| and |56| to measure the amount of entanglement in 
bipartite quantum states, one can check that the concurrence is given by 



C(t) = 2 max{0, Ai(t),A 2 (t)} 



where 



Ai(i) = 2M 2 a 2 b 2 



A 2 (t) = 2M 2 a 2 b 2 



(1 — 7)(1 + gcos??T,7r) — (1 — gcosm7r) 



(1 — 7)(1 — qcosmir) — (1 + qcosmir) 



It follows that the concurrence is given by 

at) A l ~ p2 



2 1 + p n cos mir 



e- Lt (l+p n - 2 ) - {l-p 



„n-2\ 



for 



1 



t<t = f [ ln(l+p"- 2 )-ln(l-p"- 2 )] 



and the system is entangled. However, for t > to, the concurrence is zero and the entanglement 
disappears, i.e. the system is separable. This clearly reflects that under dephasing channel, the 
entanglement suddenly vanishes. Note that the bipartite system under consideration is initially (in 
the absence of an external interaction) entangled. Indeed, for t = 0, the concurrence is 



(7(0) 



p n-2 _ p n 



1 + p n cos run 

and is always non zero except in the limiting cases p — > or p — > 1 for m even. Notice that for m odd, 
the concurrence is zero for p — > but does not vanish when p — > 1 and it is given by 2/n. 
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It is important to stress that the quantum discord Z)(pi2(0)) is nonzero except in the particular 
case p — > 0. To show that the quantum discord does not disappear after the interaction of the system 
with the dissipative channel, we note that the density matrix pi2(t)) belongs to the class of the so- 
called circulant states [37]. Thus, one can use the vanishing quantum discord criteria, discussed in [37], 
to check that the state pi2(t) has vanishing quantum discord if and only if p — > 0. This implies that 
even when entanglement suddenly disappears in a finite time, quantum discord does not vanish. This 
agrees with the commonly accepted fact that the quantum discord is more robust than entanglement 
to sudden death under a dissipative channel. 

5 Concluding remarks 

In this paper, we have obtained explicit expressions of quantum discord for symmetric and anti- 
symmetric superpositions of multipartite coherent states. These states cover the generalized GHZ 
and generalized Werner states. In particular, the balanced antisymmetric superpositions (m odd) 
interpolate continuously between generalized GHZ and generalized Werner states. The key point in 
determining the bipartite quantum correlations is based on the splitting of the entire system in two 
qubit subsystems. Two inequivalent splitting schemes were discussed. The first one leads to a pure 
bipartite density and the quantum discord coincides with the entanglement of formation. The second 
consists in constructing bipartite systems by a trace procedure keeping only the modes in which we 
are interested. In this way, mixed states are obtained and the corresponding quantum discord was 
explicitly derived. This derivation requires an optimization over all the measurement needed to extract 
the amount of quantum correlation which is general difficult to perform. To avoid such a difficulty, 
we used the purification method together with the Koashi- Winter relation which are advantageous 
in simplifying the minimization process of the conditional entropy. In the last part of the paper, we 
discussed the robustness of the quantum discord present in multipartite coherent states in comparison 
with the entanglement. We have shown that in sending the system through a dephasing channel the 
entanglement can be lost. This is not the case of quantum discord which behaves more robust against 
dissipative channels. It must be noticed that the results obtained here can be extended easily to many 
other classes of coherent states even those associated with higher symmetries which are labeled by 
several variables. Finally, we stress that the bipartite correlation does not capture genuine multipar- 
tite correlations. It follows that it is interesting to investigate the measure of genuine multipartite 
quantum discord for arbitrary multipartite non orthogonal states in the spirit of the results recently 
obtained in [58] , 
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